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FOREWORD 



The examples presented in this manual are small in size and are intended for a reader 
who is interested in an introduction to System/360 programming. While studying the 
various pieces of code, it is advisable to have access to the following manuals: 

System/360 Principles of Operation (A22-6821-1) — referred to as "OM" 

System/360 Operating Syste m — As sembly Language (C28-6514-1) — referred to as "AM" 

References are made to these manuals in the subsequent pages. 

There are three basic sections in this document. In the first two sections, the basic 
instruction set is illustrated through simple examples. It is intended that the examples 
increase in complexity, each tending to illustrate some particular point relative to a 
coding technique. Pertinent comments are included. The final section shows complete 
problems, including the necessary assembly language parameters to produce a running 
program. The problems of this section have actually been run on a machine. 

The programs shown represent "good" coding practices. They have not been optimized 
to their fullest. This optimization is left as an exercise to the reader. Where optional 
methods of approach can be taken, they are shown and discussed. 



Copies of this and other IBM publications can be obtained through IBM branch 

offices. Address comments concerning the contents of this publication to 

IBM, Technical Publications Department, 112 East Post Road, White Plains, N. Y. 10601 
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1. FIXED-POINT, LOGICAL AND BRANCH INSTRUCTIONS 

1. 1 USE OF MULTIPLE REGISTERS FOR STORING INTERMEDIATE RESULTS 

The fixed-point numbers a,b,c,d are contained in the general purpose registers GPR1, 
GPR2, GPR3, and GPR4, respectively. Calculate the quantities 

r = a+b+ | c-d | 

and s = a+b- c-d 

putting r and s in GPR5 and GPR6, respectively. 

a+b Register-to-register add. 

c-d Register-to- register subtract. 

c-d J Make the sign positive. 

r=a+b+ c-d 

- c-d J Make sign negative . 

s=a+b - c-d 



a. Fixed-point arithmetic is executed in the general purpose registers. Numbers 
in these registers are treated as signed integers (O.M. , p. 23). 

b. The above example shows the advantage of multiple accumulators. Intermediate 
results can be contained in registers, and no time-consuming reference to storage 
has to be made. (We assume that the numbers are such that no overflow occurs 
during any of the operations . ) 

1. 2 EVALUATION OF A SIMPLE ARITHMETIC EXPRESSION 

Evaluate r = -2. 0* (a+b)* d/c, where a,b, c,d are the contents of general purpose 
registers GPR1, GPR2, GPR3, and GPR4, respectively. Place the result in GPR5. 

LCR 7, 1 -a. Complement a and load into 

GPR7. 

SR 7,2 -(a+b). Subtract b from -a. 

MR 6,4 -(a+b)*d. Note multiply and divide 

conventions in MR and DR. 

DR 6,3 -(a+b)* d/c 



AR 


1,2 


SR 


3,4 


LPR 


5,3 


AR 


5,1 


LNR 


6,3 


AR 


6,1 


Comments 





LR 


5,7 


SLA 


5,1 


Comments 





Quotient to GPR5 . 
Shift to multiply by 2. 



a. This problem illustrates the positioning of operands and result in the fixed-point 
multiply and divide operations (O.M. , pp. 29-30). 

b. Multiplication or division by a power of 2 for a number in fixed-point format can 

be accomplished by shifting. This is generally more efficient than a multiplication. 
Note, in the use of shifting for division, that the 2 ! s complement notation implies a 
different rounding, as the digits shifted out are always positive. Thus a right shift 
on +25 yields +12, but a right shift on -25 yields -13. A division of -25 by 2 
yields -12, with -1 as remainder. 

c. Fixed-point arithmetic treats operands in the registers as if they were integers. 
The multiply instruction takes the multiplicand from an odd register and develops a 
double-length result in an even/odd register pair. For division the dividend is 
taken to be a double-length number situated in an even/odd register pair. The 
quotient is obtained in the odd register, while the remainder is in the even register 
(O.M. , pp. 29-30). 

1. 3 CYCLIC PERMUTATION OF A WORD GROUP USING MULTIPLE REGISTER 
COMMANDS 

Given the quantities Aj, A£. . . A^g in fullword locations starting at LOC. Cyclically 
permute the information such that it is stored in the sequence A4, A5. . ,A\q A1A2A3. 

Method 1 



LM 


0,15, LOC 


Load all registers with A^ 
through Aiq. 


STM 


0,2, LOC +52 


Store Ai, A2, A3 at end of 
LOC area. 


STM 


3, 15, LOC 


Store A4 through A^g. 


Comments 







LOC+52 refers to the byte whose address is LOC+52. Since there are four bytes per 
word, this refers to the 14th word of the group. 



Load registers with Ai,A2, A3. 
Load register with A4 through A±q. 
Store multiple in permuted order. 



Method 2 




LM 


13, 15, LOC 


LM 


0,12,LOC+12 


STM 


0,15, LOC 



Method 3 




LM 


13, 12, LOC 


STM 


0,15,LOC 


Comments 





Load registers with Ai through A±q. 
Store multiple in permuted order. 



a. Method 3 makes use of the wraparound feature of the Load Multiple instruction 
(O.M., p. 26). 

b. The above methods assume that all 16 general purpose registers can be used — 
that is, that no base registers are needed to generate address LOC. This confines 
the numerical equivalent of LOC to less than 4096. Usually, of course, some 
registers will have to be set aside as base registers, in which case the problem 
will have to be done in several steps. 

1. 4 CLEARING OF A SPACE IN STORAGE USING LOOP CONTROL TECHNIQUES 
Replace the contents of 25 fullwords starting at LOC with zeros. 
Method 1: 

Clear GPRO (see comment a). 

Clear GPR1. 

Place a 25 in GPR2. 

Clear word at LOC indexed by GPR1. 

Increment index by 4. 

Reduce count by 1 and test for completion. 

Supervisor call (see comment c). 

See comment b. 

See comment b. 



a. The fastest and most convenient way of clearing a register is to subtract it from 
itself. 

b. DC is a pseudo operation for setting up constants in the program. Hence, the steps 
COUNT and FOUR are not treated as instruction. The F denotes that the constant 
following it is expressed in 32-bit fixed-point (A.M. , p. 37). 

c. SVC is a supervisor call denoting the end of the program. This is the standard 
way of terminating the execution of a program . 





SR 


0,0 




SR 


1,1 




L 


2, COUNT 


LOOP 


ST 


0,LOC(1) 




A 


l,FOUR 




BCT 


2, LOOP 




SVC 





COUNT 


DC 


F'25' 


FOUR 


DC 


F'4' 


Comments 







Method 2 



Instead of using a BCT to control the loop, one can use a BXLE (O.M. , p. 65). 

Clear GPRO. 

Clear GPR1. 

Set GPR2 to 4. 

Set end-of-loop test to 96. 

Clear word at LOC indexed by GPR1. 

Add GPR2 to GPR1 and test whether 
greater than GPR3. 





SR 


0,0 




SR 


1,1 




L 


2, FOUR 




L 


3, END 


LOOP 


ST 


0,LOC(1) 




BXLE 


1,2, LOOP 




SVC 





FOUR 


DC 


F'4' 


END 


DC 


F'96' 


Comments 







End-of-loop test = 25x4-4. 



a. This method produces a more efficient program, but uses one more general 
purpose register. 

b. The steps L 2, FOUR and L 3, END can be replaced by LM 2,3, FOUR since the 
two words to be loaded are adjacent in memory. 



Method 3 





LM 


0,9,C 




STM 


0,9, LOC 




STM 


0,9, LOC +40 




STM 


0,4 LOC +80 




SVC 





c 


DS 


IOF'0' 


Comments 







Load multiple with zero. 

Store zeros in the first ten words. 

Store zeros in next ten words . 

Store zeros in remaining five words. 

Return to supervisor . 

Ten words of zero. 



This method is faster than the first two, but less general. 



Method 4 








SR 


0,0 




STC 


0,LOC 




MVC 


LOC(99),LOC+l 


Comments 







Set register to zero. 

Store zero into first character of field. 

Propagate zero through the 25 words 
(100 bytes). 



Method 4 is included here because of its intrinsic interest. 

1. 5 EXTRACTION OF ONE ELEMENT OF A MATRIX —TWO-DIMENSIONAL TABLE 
LOOKUP 

Given a matrix of size M x N (M rows and N columns, M,N < 1000) stored row-wise in 
consecutive fullwords beginning with A^ in MTRDC,, Given also two binary integers p 
and q in GPR1 and GPR2, respectively. Put the element Apq in GPR3. 



Load number of columns into GPR5 . 

P-l 

GPR1 contains (p-l)*N. 

q+(p-l)*N 

4*(q+(p-l)*N) Multiply by 4, using a 
left shift of 2. 

Load element into GPR3 (see comment a). 
N=5 in this example . 





L 


5,N 




S 


l,ONE 




MR 


0,5 




AR 


1,2 




SLA 


1,2 




L 


3,MTR1X-4(1) 




SVC 





N 


DC 


F'5' 


ONE 


DC 


F'l' 


Comments 







a. The element Apq has byte address MTRIX + 4*( (q-l)+(p-l)*N) or 
(MTRDC-4)+4*(q+(p-l)*N) 

1. 6 CYCLIC BIT SHIFTING 

Cyclic left shift seven bit positions the information in the 32-bit fullword LOC and 
store the result in RES. 





SR 


2,2 




L 


3,LOC 




SLDL 


2,7 




OR 


3,2 




ST 


3, RES 




SVC 





Comments 







Separate 32 bits of information into 7 
and 25 bits. 



The OR instruction is used to add the seven bits back into register 3. One could also 
use the ALR instruction to do this, but not an AR instruction. AR would give an in- 
correct result when the leading bit (sign position) of the shifted number is a 1. 

1. 7 SUM OF SQUARES USING DOUBLE-PRECISION FIXED-POINT ARITHMETIC 

Eight 32-bit integers starting at fullword boundary NUMB are given. Compute the sum 
of the squares of these integers and store it at SUM in the format sign plus a 63 -bit 
integer. It may be assumed that the sum can be expressed as a 63-bit integer. 



Set up registers . 

Load integer. 

Square the number. 

Low order of 64-bit sum . 

Branch if carry out of low-order word. 

High order of 64-bit sum . 

Close loop. 

Store result in SUM. 

Carry into high-order word by adding 1. 





LM 


0,5, REGIS 


LOOP 


L 


13,NUMB(4) 




MR 


12,13 




ALR 


3,13 




BC 


3, CARRY 


CONT 


AR 


2,12 




BXLE 


4, 0, LOOP 




STM 


2, 3, SUM 




SVC 





CARRY 


AR 


12,5 




B 


CONT 


REGIS 


DC 


F'4, 28, 0,0,0, 1' 



Comments 

a. Loading of consecutive registers can be done efficiently with a LM instruction. 

b. When a carry out of sign position in ALR instruction occurs, the condition code 
will be either 2 or 3 (O.M. , p. 27). Hence BC3 results in a successful branch 
whenever there is a carry. 

1. 8 TRANSPOSITION OF A SQUARE MATRIX USING ADDRESS CONSTANTS 

An NxM matrix has fullword elements and is stored row-wise beginning at MATRIX. 
Create the transpose of this matrix and store it in the same area. 





LM 


5, 11, PARAM 




GO 


LR 


12,9 


column pointer . 




LR 


13,10 


Row pointer . 


TRANS 


L 


3,0(12) 


Load a... 




L 


4,0(13) 


Load a... 




ST 


3,0(13) 


Store a., in a... 
i] Ji 




ST 


4,0(12) 


Store a., in a... 
Ji iJ 




AR 


12,5 


Step to next column. 




BXLE 


13, 6, TRANS 


Step to next row and cl 




AR 


9,11 


Set up for next column 




AR 


10,11 


Set up for next row. 




BCT 


8, GO 






SVC 







PARAM 


DC 


F'4' 






DC 


A(4*N,MATRK+4* 


(N*N-1)) 




DC 


A(N-1, MATRTX+4 ; 


, MATRrX+4*N,4*(N+l)) 



Comments 

a. We want to replace the element a., of the matrix by a., and vice versa. 

ij ]i 

b. DC type A is used for setting up address constants (A.M., p. 48). 

1. 9 LENGTH OF AN UNKNOWN FILE USING BYTE-HANDLING CAPABILITIES AND 
MULTIPLE REGISTERS 

Information of unknown length is written in consecutive eight-bit bytes beginning at 
INFO. Its end is signified by a special character of eight binary l's. Find the file 
length (including the special character) in bytes, and put the answer in GPR1. 



Method 1 





SR 


2,2 




L 


6, ONE 




LR 


1,6 




L 


3, SPCH 


IC 


IC 


2,INFO-l(l) 




CLR 


2,3 




BE 


EQUAL 




AR 


1,6 




BC 


15, IC 


EQUAL 


SVC 





SPCH 


DC 


X'OOOOOOFF' 


ONE 


DC 


FT 


Comments 







Set register to 0. 



GPR3 has special character. 
Load one byte of information. 
Compare it with GPR3 . 
Branch if equal (see comment b) . 

Unconditional branch (see comment b), 



a. The IC instruction can be used to handle one byte at a time. 

b. BC with mask = 15 results in an unconditional transfer. This could also be written 
in extended mnemonic form (A.M. , p. 33) simply as B. Conversely, the BE, 
which is in extended mnemonic form, could be written as BC8. 



c. DC type X is used for setting up hexadecimal constants (A.M. , p. 43). 



Method 2 





LM 


1,5, PARAM 


IC 


IC 


5,0(1,4) 




AR 


1,3 




CLR 


3,5 




BC 


7,IC 




SVC 





PARAM 


DC 


F'O' 




DC 


XL4'FF' 




DC 


F'l' 




DC 


A(INFO) 




DC 


F'O' 


Comments 







Add to counter. 

Compare with eight binary l's. 

Branch if not equal. 



a. The programmer may, if he so desires, specify which base register is to be used. 
This is shown in the second step. Alternatively, one could write the second step 
as IC 7, INFO(3) and the assembler would assign the base register. 



Method 3 










LM 


1,6, PARAM 


LOAD GP registers. 


IC 


IC 


3,0(1,6) 


Insert character into 3. 




AR 


1,4 


Increment index 




BXH 


5,2,IC 


See comments 




SVC 







PARAM 


DC 


F'O' 


Reg 1 - Index. 




DC 


F'O' 


Reg 2 - Zero. 




DC 


F'O' 


Reg 3 - Character from field 





DC 


FT 




DC 


XL4'FF' 




DC 


A(INFO) 


Comments 







Reg 4 - Index increment. 

Reg 5 - Eight l's. 

Reg 6 - Base address of field. 



The eight l's are added to zero, and then compared with the target character in 
register 3. Since eight l's will always be higher than any other character, the branch 
will be taken unless there are eight l's in the target character. Then the exit will be 
taken. This method is included to show a variation in use of index control instructions. 

Comment On A Possible Method 4 

If memory space is not critical, or if speed of function is critical, the most efficient 
and probably best method is to use the "Translate and Test" instruction. The program 
is left as an exercise for the reader. 

1.10 TABLE LOOKUP USING THE TRANSLATE INSTRUCTION 

It is desired to produce a list of eight-bit BCD code from a list of four-bit hexadecimal 
characters. (This may be used, for example, for printing.) The given list is located 
at fullword boundary HEX extending through eight doublewords for a total of 128 
hexadecimal characters. The output BCD code is to be located at fullword boundary 
BCD. 



LM 1,5, REGIS 

NEXT LH 8,HEX(1) 

LR 6,4 

LOOP SRDL 8,4 

SRL 9,4 

BCT 6, LOOP 

ST 9,BCD(5) 

AR 5,4 

BXLE 1,2, NEXT 

TR BCD(128), TABLE Translate entire line. 

SVC 



Fetch four hex characters. 
Load GPR6 with 4. 
Shift to GPR9 one hex digit. 
Expand to byte length. 

Store four expanded characters. 



10 



REGIS DC A(0, 2, 63,4,0) 

TABLE DC X'F0F1F2F3F4F5F6F7' 

DC X'F8F9C1C2C3C4C5C6' 

Comments 

The translate instruction (TR) works like a table-lookup scheme. It translates 256 bytes 
with just one setup. Here the eight-bit arguments are located at BCD and the eight-bit 
function bytes are located at TABLE. The translation of each argument byte is com- 
pleted with the function byte replacing the argument at BCD. 
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2. FLOATING-POINT INSTRUCTIONS 
2. 1 POLYNOMIAL EVALUATION 

Evaluate the polynomial 

20 
P(x) = y* Akx k 

K=0 

where x is located at X, a 2 Q at A, a^9 at A+4, etc. Store the result (single precision) 
in POLY, 





SR 


1,1 




LM 


2, 3, INCRE 




LE 


o,x 




SER 


2,2 


LOOP 


AE 


2,A(1) 




MER 


2,0 




BXLE 


1,2, LOOP 




AE 


2,A(1) 




STE 


2, POLY 




SVC 





INCRE 


DC 


F'4' 




DC 


F'76' 


Comments 







Clear GPR1. 

Load index value and limit. 

Load argument in FPRO. 

Zero FP register 2. 

Add coefficient of kth term. 

Multiply by x. 

Add coefficient of 0th order. 

Store result. 

Return to supervisor. 



The above technique is known as "nesting". P(x) =(.... (A 2 qX + A 19 ) x +. 
technique is more efficient than a term by term evaluation. 



+ A-, ) X + Ar 



2. 2 SEPARATION OF INTEGER AND FRACTION PARTS OF FLOATING-POINT 
NUMBERS USING UNNORMALIZED INSTRUCTIONS 

The single-precision (32-bit) floating-point number N in location DOG has a small 
(4 6) exponent magnitude. Create two floating-point numbers I, F in CAT and CAT+4 
such that 



and 



I = an integer 
| F| < 1. 0, sign of F = sign of N 

1+ F = N. 



SDR 


6,6 


LDR 


2,6 


LE 


6, DOG 


AW 


6,X6 


LER 


2,6 


SDR 


6,2 


AU 


6,X0 


STE 


2, CAT 


STE 


6, CAT+4 



Zero FPR 6 (double). 

Zero FPR 2 (double). 

Load N. 

Integer in 1st and fraction 2nd half 

of FPR6. 
Load integer into FPR2. 
Fraction in FPR6. 
To force an unnormalized fraction. 
Store integer part. 
Store fractional part. 
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SVC 







DC 


OD 


X6 


DC 


X'4600000000000000' 




DC 


X'40000000' 


Comments 







Place X6 on doubleword boundary. 
Exp 6, fraction 0. 
Exp 0, fraction 0. 



a. I and F are generally unnormalized, the exponent of I being 6 and that of F zero. 

b. Unnormalized addition of a number with a zero fraction can be used to force the 
exponent of number to a predetermined value. 

2. 3 EXTRACTING THE INTEGER PART OF A FLOATING-POINT NUMBER 

The 32-bit single-length floating-point number N in location DOG has a small exponent 
( <6) and is therefore less than 224 i n magnitude. Put the integer part of the number 
in GPR1 in fixed-point integer form. 

Method 1 



Integer part with exponent of 6 
in FPR4. 

Integer part with exponent in GPR1, 
Blank out exp. , leaving sign. 
Test whether number is negative. 
Branch if positive or zero. 
Remove sign. 
Take 2's complement. 



Mask to eliminate exponent. 
Floating zero with 6 exponent. 
Mask to eliminate sign. 



A negative number in fixed point is represented in two's complement form, while in 
floating-point the fraction is represented in true form. Hence special steps have to be 
taken to create the complement if necessary. 

Method 2 





LE 


4, DOG 




AU 


4,X6 




STE 


4, TEMP 




L 


1,TEMP 




N 


1, XZERO 




LTR 


1,1 




BC 


10, NOTNEG 




N 


1, NSIGN 




LNR 


1,1 


NOTNEG 


SVC 





TEMP 


DS 


IE 


XZERO 


DC 


X'80FFFFFF' 


X6 


DC 


X'46000000' 


NSIGN 


DC 


X'7FFFFFFF» 


Comments 







This method does not 


use the fixed- 


-point registers, 




LE 


4, DOG 




Load N into FPR 4. 


AU 


4,X6 




Integer part with exponent of 6 in 
FPR4. 


BC 


10, NOTNEG 


Branch if positive or zero. 


AU 


4, ONES 




Form l's complement. 


AU 


4, ONE 




Form 2's complement. 


STE 


4,1 




Store result in I. 


MVI 


I,X'FF' 




Place FF in exponent. 


SVC 









NOTNEG STE 


4,1 




Store result in I. 
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MVI 
SVC 


1,X'00' 



I 

X6 
ONES 


DS 
DC 
DC 


F 

X'46000000' 

X T 46FFFFFF' 


ONE 


DC 


X'46000001' 


Method 3 






EXIT 

TEMP 

MASK 


LE 

AW 

STD 

L 

BC 

LCR 

SVC 

DS 

DC 


4, DOG 

4, MASK 

4, TEMP 

l,TEMP+4 

10, EXIT 

1,1 



ID 

X'4E00000000000000» 


Comments 







Blank out exponent. 

Align on fullword boundary 

Unnormalized zero. 

Maximum positive floating-point 

integer. 
Unnormalized 1. 



Integer part to right end of FPR4. 
Does not change condition code. 
Does not change condition code. 



a. The AW instruction forces the integer part of the number to the extreme right end of 
the 64-bit floating-point register. 

b. The condition code set by the AW instruction is not affected by the two subsequent 
instructions. 

Method 4 



LE 
AU 

STE 

L 

LPR 



SLDA 



4, DOG 
4,X6 

4, TEMP 
0,TEMP 
1,0 



0,39 





SRA 


0,7 




SVC 





TEMP 


DS 


IE 


X6 


DC 


X»46000000» 


Comments 







Integer part with exponent of 6 in 
FPR4. 

Integer part with exponent in GPRO. 
Load absolute value if positive and 

complement if negative into GPR1. 

Correct sign still in GPRO. 
Shift quantity against sign and shift 

out the exponent. 
Shift to proper position. 



LPR creates the correct quantity except for the exponent part, now extraneous, and sign. 
The left double shift removes the exponent and regains the sign. The right shift provides 
the proper offset. Negative numbers will have leading 1-bits as prescribed by the two's 
complement notation. 
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2. 4 CONVERTING A FIXED-POINT NUMBER TO FLOATING-POINT 

GPRO contains a small integer which is less than 2^ in magnitude. The integer is to be 
converted into a floating-point number and left in FPRO in double-precision form. 



Method 1 





LR 


4,0 




N 


4, SMASK 




LPR 


3,0 




OR 


3,4 




O 


3,X6 




ST 


3, TEMP 




LD 


0,TEMP 




AE 


0,X6 




SVC 





TEMP 


DS 


ID 


SMASK 


DC 


X' 80000000' 


X6 


DC 


X'46000000' 


Comments 







Separate sign bit and leave in GPR4, 

Load absolute value into GPR3. 

Put back sign bit. 

Tag on exponent of 6. 

Store into left half of TEMP. 

Load doubleword with right half of 

FPRO zero. 
Normalize the number. 



The right half of FPRO may have contained extraneous bits, 
doubleword avoids an extra zeroing instruction. 



The use of TEMP as a 



Method 2 





LPR 


1,0 




O 


1,MASK 




SLDA 


0,32 




ST 


0,TEMP 




LD 


0,TEMP 




AE 


0,MASK 




SVC 





TEMP 


DC 


D'0. 0' 


MASK 


DC 


X46000000 


Comments 







Load absolute value into GPR1. 
Insert an exponent of 6. 
Shift fraction and exponent into 
GPRO. Original sign is preserved. 



To force normalization. 



a. The LPR instruction creates the absolute value (with eight lead zero-bits) in GPR1, 
the SLDA instruction shifts the GPR(0, 1) pair to result in the correct sign-magnitude 
notation. The additional instructions move the result to FPRO and perform normal- 
ization. 

b. AE can be used to force normalization since the right-hand side of FPRO is known to 
be zero. 
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2. 5 DOUBLE-PRECISION SQUARE ROOT SUBROUTINE 



Given a double-precision number x in FPRO, create the double-precision square root 

in FPRO 

Algorithm: let x=16 . f wherej^- < f<l 

then an approximation vq to y = Jx can be obtained by the linear expression 

yo = (*L + l\ [E/2] 

■ y ° \9 9/16 

/32f 8\ [E/2] 

y o = (,— + 7 j 16 

(Here [e/2] denotes the integer part of E/2.) 

It can be shown that such an approximation has a relative error which is nowhere greater 
than 1/9. Therefore, four Newton-Raphson iterations of the form 



when I is even 



when I is odd 



* ( V V 



n+ 1 2 

are sufficient to give better than 56-bit accuracy. 

Conventions 

In order to gain speed, the loop for the Newton-Raphson iterations was expanded and 
single precision used whenever possible. In order to achieve the multiplication by 1/2, 
the HALVE instruction is used when possible. However, this cannot be done on the 
last step, since HALVE produces an unnormalized result. This can be overcome by 
using a multiply by 1/2 on the last step. Alternatively one could use 



2 \ Y 3 s ) 



Y„ = Y + — 
4 3 



which is faster than a multiply but may give less accuracy (extra roundoff). The 
argument is to be given in FPRO and the answer is returned in FPRO. GPRO, GPR1, 
FPRO, FPR2, and FPR4 are used. 



SQRTD 



LTDR 


2,0 


BC 


12, EXCP 


STE 


0,TEMP 


SR 


0,0 


IC 


0,TEMP 


SRL 


0,1 


A 


0,MSKS1 


SR 


1,1 


TM 


TEMP, X' 01' 


BC 


8, EVEN 



Test argument. 
Arg. or negative. 

Zero GPRO. 

Insert exponent byte in GPRO. 

Divide exponent by 2. 
E/2 Regenerate excess -64 charac- 
teristic. 

Zero GPR1. 

Test exponent for even/odd. 

Branch to even if exponent even. 
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ODD 


L 


l,FOUR 


If exponent odd, set index to 4. 


EVEN 


STC 


0,TEMP 


Insert [e/2] 




STC 


0,C2(1) 


Insert [e/2J 




LE 


0,TEMP 


Load f [E/2J 




ME 


0,C1(1) 


Multiply by either 8/9 or 32/9. 




AE 


0,C2(1) 


Initial approximation. 




LER 


4,2 


Start of 1st iteration. 




DER 


4,0 


x/Y n 




AER 


0,4 


(Y n + x/Y n ) 




HER 


0,0 


l/2(Y n + x/Y n ) 




LER 


4,2 


Start of 2nd iteration. 




DER 


4,0 






AER 


0,4 






HER 


0,0 






LDR 


4,2 


Start of 3rd iteration-double pre 
cision. 




DDR 


4,0 






ADR 


0,4 






HDR 


0,0 






DDR 


2,0 


Start of 4th iteration. 




ADR 


0,2 






MD 


0,HALF 


Final result. 




B 


0(14) 


Normal return. 


EXCP 


BC 


4, ERROR 


Error return. 




B 


0(14) 


Normal return for zero arg. 


TEMP 


DS 


IE 




MSKS1 


DC 


X»00000020' 




FOUR 


DC 


X , 00000004« 




C2 


DC 


E'.222» 


2/9 




DC 


E'.8889' 


8/9 


CI 


DC 


£'.8889' 


8/9 




DC 


E T 3. 55556' 


32/9 


HALF 


DC 


D'.5» 




Comments 









DC type D is used for setting up long floating-point constants. The constant is aligned 
at the proper doubleword boundary (A.M., p. 45). 

2.6 MATRIX ADDITION 

Two NxM matrices, whose elements are single-length (32-bit) floating-point numbers 
stored row-wise starting at AMTX and BMTX respectively, are to be added and the 
resultant matrix is to be stored row -wise starting at RMTX. 
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Method 1 



LM 
SR 



SR 



1,5, REGIS 
1,3 



2,3 



Load addresses. 

Form difference so that when added 

to GPR3 the correct AMTX element 

address is formed. 
Form difference for BMTX (see 

above) . 



LOOP 


LE 


0,0(1,3) 


Load AMTX element. 




AE 


0,0(2,3) 


Add BMTX element. 




STE 


0,0(3) 


Store in RMTX. 




BXLE 


3, 4, LOOP 






SVC 







REGIS 


DC 


A(AMTX, BMTX, RMTX, 


4, RMTX+4(M*N-1) 


N 


EQU 


3 


N=3 for this example. 


M 


EQU 


5 


M=5 for this example. 


Comments 









a. Requires three registers for indexing but makes no restrictions on N and M. 

b. EQU is used to define a symbol (A.M., p. 36). 

c. RMTX+4*(M*N-1) is the address of the last element on RMTX. 
Method 2 



Address of AMTX to GPR1 
Address of BMTX to GPR2 
Address of RMTX to GPR3 





L 


5, LIMIT 




SR 


13,13 




LA 


1,AMTX 




LA 


2, BMTX 




LA 


3, RMTX 


LOOP 


LE 


0,0(13,1) 




AE 


0,0(13,2) 




STE 


0,0(13,3) 




BXLE 


13, 4, LOOP 




SVC 





LIMIT 


DC 


A(4*(M*N-1)) 



where N and M would be specified and GPR13 would take on values to 4* (M*N-1) in 
increments of 4. 

Comments 

a. Requires four registers for indexing. 
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Method 3 






LOOP 


LE 


0,AMTX(13) 




AE 


0, BMTX(13) 




STE 


0, RMTX(13) 




BXLE 


13, 4, LOOP 


Comments 







a. Uses only one register for indexing but may not work when the matrices are large. 

b. Another method would be to intersperse the matrices. 

2. 7 PRODUCT OF TWO SQUARE MATRICES 

Two NXN single-precision (32-bit) floating-point matrices L and R are stored row-wise 
beginning at LMTX and RMTX respectively. Create P=L*R and store it row-wise 
beginning at PMTX. ( Note: C(12) means the contents of register 12.) 



* 
* 
* 
LOOP 2 



LOOP 1 



N 

ROW 
SQMX 
REGIS 



LM 



LR 

LR 

SDR 

LE 

ME 

AER 

AR 

BXLE 

STE 

AR 

BXLE 

SR 

AR 

BXLE 

SVC 

EQU 

EQU 

EQU 

DC 

DC 



6, 14, REGIS 



3,13 

2,12 

0,0 

2,0(2) 

2,0(3) 

0,2 

3,8 

2, 6 LOOP1 

0,0(14) 

14,6 

13,10,LOOP2 

13,8 

7,8 

12, 8, LOOP2 



1000 

4*(N-1) 

4*(N*N-1) Row +4 

A(4, LMTX + ROW, ROW + 4 LMTX + SQMX) 

A(4, RMTX+ROW, LMTX, RMTX, PMTX) 



C(12) = Starting address of a row of 

LMTX. 
C(13) = Address of RMTX. 
C(3) = Index down columns of RMTX. 
C(2) = Index across rows of LMTX. 
C(FPRO) = Sum of products. 
Start of inner loop. 
Form product of elements. 
Running sum. 
Step to next row. 
C(2) = C(2)+4 and loop. 
Store element. 
C(14) = C(14)+4. 
C(13) = C(13)+4 and loop. 
C(13) = address of RMTX. 
Move to next row of LMTX. 
C(12)=C(12)+4*N and loop. 

Any positive integer. 



Comments 

a. An asterisk in column 1 indicates a comments card (A. M. , p. 11). 
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3. COMPLETE PROBLEMS 



3. 1 SORTING ON THE BASIS OF SUBFIELDS (see comments) 



Given 256 consecutive fields each having an "A" subfield of eight bits and a "B" 
subfield of 16 bits, sort on the basis of "A" at ANS. (All "A" fields are different. ) 





A 


B 






(8 bits) 


(16 bits) 






ORG 


4096 




BEGIN 


BALR 


1,0 






USING 


M 






LM 


4, 7, REGIS 




LOOP 


LM 


9,11,0(6) 






SR 


12,12 






IC 


12,9(,6) 






AR 


12,12 






STH 


11,0(12,7) 






SRDL 


10,24 






SR 


12,12 






IC 


12,6(,6) 






AR 


12,12 






STH 


11,0(12,7) 






SRDL 


10,24 






SR 


12,12 






IC 


12,3(,6) 






AR 


12,12 






STH 


11,0(12,7) 






SRL 


9,8 






SR 


12,12 






IC 


12,0(,6) 






AR 


12,12 






STH 


9,0(12,7) 






BXLE 


6,4, LOOP 






SVC 







ANS 


DS 


OF 






DS 


256H 




REGIS 


DC 


A (12 , DAT A+756 , DAT^ 


l,ANS) 


DATA 


DC 


X' 

256 data fields 






END, BEGIN 





See comment a. 

See comment b. 

Process four fields at a time. 

Clear GPR12. 

Fetch "A" field into GPR12 

(note displacement of nine 

bytes). 

Double contents of GPR12. 

Store "B" field (two bytes). 

Shift to get B at any even boundary 



Store second field. 



Store third field. 



Store fourth field. 

Load to process four more 

fields . 

See comment d. 



See comment c. 
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Comments 

a. ORG is a pseudo instruction defining the beginning address of the program. 

b. Assignment of base registers is done automatically by the assembler, but the 
programmer must: 

1. Specify what registers may be used as base registers and inform the assembler 
of their contents 

2. Load the base registers with the appropriate values 
The standard method of achieving this is through the sequence 



BALR 
USING*, R 



R, 

*,R 



Where R is a register. 



c. END BEGIN is a pseudo instruction defining the end of the program. It also tells 
the monitor to start execution at BEGIN when assembly is completed. 

d. Word alignment on the fullword boundary can be achieved through a DS OF pseudo 
instruction (A.M. , p. 51). 

e. This problem shows how information whose field length is not a fullword multiple 
can be handled. The reader is urged to work through this problem, showing the 
contents of the register at each step. 

3. 2 REMOVAL OF KEYWORDS 



Given a string of 100 eight -bit bytes at DATA to DATA+99 (25 fullwords). Remove the 
words which match the four -letter keyword KEY and place the condensed result at ANS. 

Method 1: Character-by-Character Operation 





ORG 


4096 


START 


BALR 


15,0 




USING 


*,15 




LA 


3 , DATA 




LA 


4, ANS 




L 


0,ONE 




L 


1,NINE6 




AR 


1,3 




L 


6, KEY 




L 


9, DATA 




B 


CMP 


EQUAL 


A 


3, FOUR 




CR 


3,1 




BH 


HIGH 


LODE 


IC 


9,0(,3) 



Address of data to GPR3. 
Address of answer to GPR4. 



Ending address of data field. 
First four characters. 



Exit if data exhausted. 
Load one byte. 
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SLL 



9,8 





IC 


9,1(,3) 




SLL 


9,8 




IC 


9,2(,3) 




SLL 


9,8 




IC 


9,3(,3) 




B 


CMP 


IC 


IC 


9,3(,3) 


CMP 


CLR 


6,9 




BE 


EQUAL 




SLDL 


8,8 




STC 


8,0(,4) 




AR 


4,0 




BXLE 


3,0,IC 


HIGH 


A 


l,FOUR 




SR 


1,3 




BZ 


END 




STC 


1,MVC+1 


MVC 


MVC 


0,(1,4), 0(3) 


END 


SVC 





ONE 


DC 


FT 


FOUR 


DC 


F'4' 


NINE 6 


DC 


F'96' 


ANSW 


DS 


25F 


KEY 


DC 


C'ABC » 


DATA 


DC 


C'ABC this ABC is 




END 


START 


Comments 







Shift eight positions to make 

room. 

Load another byte. 

Shift eight positions. 

Load 

Shift 

Load 

Load another byte. 

Compare four bytes to KEY. 

Branch if equal. 

Shift out one byte. 

Store the byte. 

Step answer field pointer by 1. 

Increment GPR3 by 1. 



Move the last bytes. 



Keyword (note that last 
character is a blank). 



The MVC instruction has an SS format. The assembler language format for these 
instructions differs from previously encountered format (A.M. , p. 24). 

The instructions between LODE and the subsequent branch could have been replaced by 



LODE 



MVC 


LOC,0(3) 


L 


9,LOC 


B 


CMP 



LOC 



DC 



F'O' 
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Method 2 

Use the TRT instruction to detect the occurrence of leading character in KEY, then 
examine next three characters for matchedness. 





ORG 


4096 




START 


BALR 


15,0 






USING 


*,15 


Rl START ADDR &END OF HIT 
(START- 1). 
INITIALIZE R2 FUNCT=2. 




LM 


-< 
1,6, HOSKP 


WORKING R3 START+98. 

REGISTERS R4 0, NEW START ADDR WORK REG. 
R5 0, LENGTH OF REMAINING SCAN 
L . R6 ANSWER ADDRESS. 


LOOP2 


LR 


4,1 


New start addr— *-R4. 


LOOP1 


LR 


5,3 


End addr— >»-R5. 




SR 


5,1 


Length of TRT. 




EX 


5, TRT 


Execute TRT inst. 




BC 


5, OUT 


End of field. 




CLC 


KEY+1(3), 1(1) 


Test remainder of key. 




BNE 


LOOP1 


Continue test cycle on ^. 




SR 
SR 


1,4 
1,2 


Set up length for data move to answer area. 




EX 


1,MVC 


Execute data move. 




LA 
LA 


6,1(1,6) 
1,5(1,4) 


, Set up to test remainder of the field. 




BC 


LOOP2 


Continue test. 


OUT 


SR 


3,4 


Calc. final length. 




EX 


3, MVC 


Execute final move. 




SVC 





Return to monitor . 


TRT 


TRT 


1(0, 1), TABLE 




MVC 


MVC 


0(0,6), 1(4) 




HOSKP 


DC 


(Storage for 6 


reg. loading) 


TABLE 


DC 


(Table for TRT inst (Function=2)) 



PRIME NUMBER GENERATOR 

Find all prime numbers whose value is less than a certain number LIM. Algori thm : To 
decide whether a number N is a prime, we divide it by all primes with value ^: ■yN. If 
no exact divisor is found, N is a prime. 



START 



ORG 


4096 






BALR 
USING 


15,0 
*,15 






L 


14, ONE 






ST 


14, M 


First number = 


= 1 


AR 


14,14 






ST 


14,M+4 


Second prime = 


= 2 


A 


14, ONE 
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ST 


14,M+8 




SR 


13,13 




ST 


13, MNS 




ST 


14, NM 




L 


9, TWELVE 


LOOP1 


L 


14, NM 




A 


14.M+4 




ST 


14, NM 




S 


14, LIM 




BP 


STOP 




L 


13 , EIGHT 


LOOP2 


LM 


10, 11, MNS 




D 


10,M(13) 




BXLE 


10,11,LOOP1 




S 


11,M(13) 




BM 


PRIME 




A 


13, FOUR 




B 


LOOP2 


PRIME 


L 


8,NM 




ST 


8,M(9) 




A 


9, FOUR 




B 


LOOP1 


STOP 


SVC 





MNS 


DS 


IF 


NM 


DS 


IF 


ONE 


DC 


X T 00000001' 


FOUR 


DC 


X'00000004' 


EIGHT 


DC 


X' 00000008' 


TWELVE 


DC 


X'OOOOOOOC 


LIM 


DC 


X' 00000064' 


M 


DS 


50E 




END 


START 



Third prime = 3. 

Zero GPR13. 

Zero MNS 

Store C(14) in previous prime number. 

GPR9 = Index on prime number array. 

Load previous number. 

Add 2 to previous number. 

Next number to be tested. 

Test for end. 



Load current numbers to be tested. 

Remainder in Reg. 10. 

If Rem = 0, number is not prime, so branch to 

LOOP1 to next candidate. 

Square root test. 

Number is prime. 

Try next prime number as divisor. 

Load number found to be prime. 
Store in prime array. 
Bump prime array index. 



Previous prime number. 



LIM = 100 in this example. 
Prime numbers stored here. 



3.4 MIN- MAX PROBLEM 

Given an array of N rows and M columns, whose elements are 32-bit integers stored 
row -wise beginning at ARRAY. Find the maximum of each row and store the addresses 
of these maxima in fullwords beginning at MAX. Then find the minimum of these 
maxima and store its address at MIN. For this example an array of three rows and 
five columns is used. 



ORG 



4096 



START 


BALR 


15,0 




USING 


*,15 




LM 


2, 7, REGIS 




L 


1,0(3) 




ST 


3, MIN 


REPLACE 


L 


14,0(3) 



C(l) = MIN of MAX. 

Address of current minimum. 

C(14) = MAX. 
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ST 


3,MAX(2) 




Address of current maximum for this row. 




B 


BXH 






CMP 


C 


14,0.(3) 




Compare current maximum with operand nev 




BL 


REPLACE 




maximum discovered. 


BXH 


BXH 


3, 4, CMP 




Repeat until row is exhausted. 




CR 


1,14 




Is NEW MAX less than MIN? 




BH 


NEWMIN 




Branch if so. 


AR 


AR 


2,4 




Handle next member of MAX vector. 




BXH 


5, 6, REPLACE 


Proceed to next row. 




SVC 









NEWMIN 


LR 
L 


1,14 
10,MAX(2) 




New minimum generated. 




ST 


10, MIN 




Address of new MIN saved. 




B 


AR 






N 


EQU 


3 






M 


EQU 


5 






END 


EQU 


4*(N*M-1) 






MT4 


EQU 


M*4 






ARRAY 


DC 


F'1,3,5,7,. 


-9,-1, 


-3,-50,-7,8,2,4,0,8,10' 


MAX 


DS 


5E 






MIN 


DS 


IE 






REGIS 


DC 


A(4*N-4, ARRAY+END , -4, ARRAY+END-MT4, -MT4, ARRAY-8) 




END 


START 






Comments 











a. Candidate for max is at GPR14, with address of candidate at MAX(2). Candidate 
for min is at 1, with address at MIN. 

3.5 DOUBLE-PRECISION BINARY- TO-DECIMAL CONVERSION 

Given a 64-bit binary positive integer beginning at fullword boundary BIN, convert it 
into decimal beginning at DEC. 



START 



CONT 
ADD 



ORG 


4096 


BALR 


15,0 


USING 


*,15 


LM 


2, 3, BIN 


CVD 


2,DEC+8 


MP 


DEC(16),TW032(6) 


SLDA 


2,31 


CVD 


2, TEMP 


MP 


TEMP(8),TWO(l) 


LTR 


3,3 


BC 


4, ADD 


AP 


DEC(16),TEMP(8) 


SVC 





AP 


TEMP(8),ONE(l) 


B 


CONT 



64-bit number to GPR2 and GPR3. 

Convert high-order part. 

Decimal-multiply high-order by 2 . 

Shift low order to GPR2. Save low -order bit 

in sign of GPR3. 

Convert low-order part. 

Decimal-multiply low order by 2. 

Set condition code. 

Test for low-order bit. 

Add two parts together. 

Add in low -order one. 
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TW032 


DC 


X' 0429496 7296C 


TWO 


DC 


X'2C 


ONE 


DC 


X'lC 


TEMP 


DS 


ID 


DEC 


DS 


2D 




END 


START 


Comments 







Decimal constant 2 3 ^. 
Decimal constant 2-^. 
Decimal constant 1. 



a. A 63-bit binary integer plus sign is equivalent to a decimal number of not more 
than 20 digits plus sign. A 128-bit field is used to store this decimal result. 

b. The CVD instruction converts a 32-bit signed integer to decimal. The 63-bit 
magnitude field is, therefore, broken into two 31-bit fields and a 1-bit field. The 
high-order part is converted and then multiplied by 2^2; the lower part is converted 
and multiplied by 2. The remaining bit contributes either a or 1 to the result. 

3. 6 SCANNING A MESSAGE FOR SYNTACTICAL ERRORS 

SIMSCAN is a very simple scanner which could be used for precompilation scanning of 
programs. In this use it would locate low-level syntactic errors without wasting the 
time of the compiler proper. 

The code shows a standard use of TRT (translate and test). It also exemplifies use of 
several general registers for keeping track of a multidimensional situation. 

Problem 

Given a string of characters starting at a known location and terminating with a period. 
The string is called error-free if and only if: 

1. There is either one or no equal sign. 

2. No two arithmetic connectives (+,-,*,/) are adjacent. 

3. At no point during a left-to-right scan have more right parentheses than left 
parentheses been encountered, and the total number of occurrences of each is the 
same. 

Solution 

SIMSCAN is written as a quasi-subroutine. The program below contains not only 
SIMSCAN but a main program called CONTROL and several input strings. 
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The register usage table, which follows, explains the code: 



General Purpose 
Register No. 



Contents 



CONTROL 



NEWDATA 



FINISHED 
ERROR 



ADSCAN 

ADSTRNGS 

FOUR 





1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 

CSECT 

BALR 

USING 

L 

L 

L 

L 

L 

L 

SR 

AR 

LR 

BR 

AR 

L 

LTR 

BM 

AR 

LR 

BR 

SVC 

SVC 

DC 

B 

DC 

DC 

DC 



= 1 

Argument address for TRT, 
Function byte for TRT. 
(not used) 



= 256 

Left -right parenthesis count (start at 0). 

Equal sign count (start at -1). 

Address of last previous connective encountered. 

Base address for STRINGS. 

Base address for SIMSCAN. 

= 4 

Length of current statement. 

Program base address. 



15,0 
*,15 
7,TW056 

12, ADSCAN 
11, ADSTRNGS 
14,0(11) 

13, FOUR 
0,ONE 
2,2 
11,13 
1,11 

12 

11,14 

14,0(11) 

14,14 

FINISHED 

11,13 

1,11 

12 



1 

X'0010' 

NEWDATA 

A(SCAN) 

A(STRINGS) 

F'4' 



Controlling program. 
Addressability. 



Load bases. 

Length of current string. 
= 4 
= 1 

Skip length word. 
Input to simscan. 
Call simscan. 
Address of new input. 

Check for end signal 



Stop 
Print 
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ONE 


DC 


F'l' 




TW056 


DC 


F'256' 




SIMSCAN 


CSECT 




* * * Scanner * * * 




USING 


*,12 


Loaded by control. 




SR 


8,8 


Clear scratch pad. 




LCR 


9,0 


Clear scratch pad. 




LCR 


10,13 


Clear scratch pad. 




SR 


1.0 


Initialize. 


CONTINUE 


AR 


1.0 


Advance beyond detected symbol 




TRT 


0(255.1), TABLE 


Start or continue scan. 




BC 


6,SIMSCAN(2) 


Symbol found, take action. 




AR 


1,7 


No symbol, advance, 




B 


CONTINUE 


Continue. 


LPAREN 


AR 


8,0 


Increase (count. 




B 


CONTINUE 




RPAREN 


SR 


8,0 


Decrease (count. 




BM 


ERROR 






B 


CONTINUE 




EQUALS 


AR 


9,0 


Increase = count. 




BP 


ERROR 


Check for . GT. 1. 




B 


CONTINUE 




CONNECT 


AR 


10,0 


See if last previous 




SR 


10,1 


connective 




BC 


10, ERROR 


was too close. 




LR 


10,1 


Reset 




B 


CONTINUE 




PERIOD 


LTR 


8,8 


End of message 




BP 


ERROR 


Check (count. 




B 


NEWDATA 




TABLE 


DC 


9D'0' 


for ASCII code 




DC 


ALl(LPAREN-SIMSCAN) 






DC 


ALl(RPAREN-SIMSCAN) 






DC 


ALl(CONNECT -SIMSCAN) 






DC 


AL1 (CONNECT -SIMSCAN) 






DC 


2X'0' 






DC 


ALl(CONNECT -SIMSCAN) 






DC 


ALl(PERIOD-SIMSCAN) 






DC 


ALl(CONNECT-SIMSCAN) 






DC 


13X'0' 






DC 


AL1 (EQUALS-SIM SCAN) 






DC 


X»00' 






DC 


20D'0' 




STRINGS 


CSECT 








DC 


F'30' 






DC 


C'THIS (E*F-L+-) 


i 




DC 


F'31' 






DC 


C'BEEP (())()(() )) )(( 


t 




DC 


F'10' 






DC 


C'X=Y+Z=2 ' 






DC 


F'10' 





28 





DC 


C'((X)((. 




DC 


F'10' 




DC 


C'X+Y(=)-2Z.» 




DC 


F'-l' 


TABLE 


DC 


9D'0 f 




DC 


CL3 T 000' 




DC 


ALl(PERIOD-SIMSCAN) 




DC 


CLl'O' 




DC 


ALl(LPAREN-SIMSCAN) 




DC 


ALl(CONNECT-SIMSCAN) 




DC 


18CL1'0' 




DC 


ALl(CONNECT-SIMSCAN) 




DC 


12CL1'0' 




DC 


ALl(CONNECT-SIMSCAN) 




DC 


15CL1'0' 




DC 


AL1 (EQUALS-SIM SCAN) 




DC 


CLl'O' 




DC 


16D'0' 



for EBCDIC code 



29 



C20- 1620-0 



S3 
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